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Motivated by potential realizations in cold-atom or cold-molecule systems, we have performed 
quantum Monte Carlo simulations of two-component gases of fermions in two dimensions with hard- 
core interactions. We have determined the gross features of the zero-temperature phase diagram, by 
investigating the relative stabilities of paramagnetic and ferromagnetic fluids and crystals. We have 
also examined the effect of including a pairwise, long-range r~^ potential between the particles. Our 
most important conclusion is that there is no region of stability for a ferromagnetic fluid phase, even 
if the long-range interaction is present. We also present results for the pair-correlation function, 
static structure factor, and momentum density of two-dimensional hard-core fluids. 

PACS numbers: 67.85.Lm, 02.70.Ss, 67.10.Db 
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I. INTRODUCTION 

Ultracold gases of fermionic atoms and molecules have 
been the subject of a large number of experimental and 
theoretical studies^i^ in the last decade. The interest 
in these systems arises from the fact that the atoms 
or molecules obey quantum statistics rather than clas- 
sical statistics and exhibit interesting and novel quan- 
tum phases at sufficiently low temperatures. Unlike elec- 
tron systems, it is possible to manipulate the interac- 
tion between the atoms or molecules to some extent, 
e.g., by Feshbach resonances^ and/or by applying mi- 
crowave fields 4r— The ability to control the interactions 
can be used to gain insight into the physics of phase tran- 
sitions and the nature of correlation effects in these quan- 
tum phases. Celebrated examples are the experimental 
studies of the Mott transition for bosons^i and of the su- 
perfluid pairing in the unitary (strong-coupling) Fermi 
gasi^ Finally, ultracold atomic gases may play an im- 
portant role in future quantum computing devices. 

Recently, interest has turned towards the use of ul- 
tracold atomic systems to investigate ferromagnetism in 
Fermi gases. Itinerant ferromagnetism in electron sys- 
tems is poorly understood, and it is possible that in- 
sights might be gained by studying ferromagnetic fluid 
states in cold-atom systemsi^ Experimental studies of 
strongly repulsive Fermi gases close to Feshbach reso- 
nances have found behavior consistent with the forma- 
tion of (local) ferromagnetic correlations in the manner 



of Stoner ferromagnetism.— However, these systems are 
limited by an intrinsic instability to molecule formation, 
which results in a short lifetime for the ferromagnetic 
state and may also complicate the interpretation of ex- 
perimental features . ^^'^^ It is therefore of interest to ex- 
plore other forms of interaction that might give rise to 
ferromagnetism in the absence of nearby bound states. 

An interesting class of experimental system in this re- 
gard is provided by gases of fermionic polar molecules 
(e.g., ^^K^'^Rh or "^W^^K). Polar molecules that are con- 
fined to two-dimensional (2D) layers and dressed by a 
circularly polarized microwave field experience a long- 
range potential that falls off as r~^, where r is the inter- 
molecular distancei^"— Furthermore, the dressed states 
have the feature that there is a very sharp increase in 
the potential energy at short distances, preventing close 
approach. Hence such molecules can be modeled as a 
gas of fermionic hard-core particles with an additional 
long-range potential varying as . The hard-core po- 
tential prevents the formation of shallow bound states 
even for (weakly) attractive long-range interactions. It 
has been suggested that a single-component gas of such 
fermionic molecules will exhibit a topological superfluid 
p-wave pairing phase, of direct relevance to possible topo- 
logically protected quantum computing devices.— 

Here we consider the case of a two-component gas 
of such fermionic molecules. We have applied quan- 
tum Monte Carlo (QMC) methods to address the issue 
of whether or not the hard-core interactions cause this 
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system to exhibit a region of itinerant ferromagnetism. 
In particular we have used the variational and diffusion 
quantum Monte Carlo (VMC and DMC) techniquesii^i^ 
to establish the ground-state phase diagram of two- 
component gases of fermionic hard-core particles mov- 
ing in 2D. QMC methods are widely acknowledged to be 
the most accurate first-principles techniques available for 
studying condensed matter. In the VMC method, Monte 
Carlo integration is used to evaluate expectation values 
with respect to an explicitly correlated many-body trial 
wave function. In the DMC method, the ground-state 
component of a trial wave function is projected out by a 
stochastic algorithm. Fermionic symmetry is maintained 
by the fixed-node approximation.— 

Some of the earliest applications of QMC methods were 
to study the ground-state properties of bosonic hard- 
sphere fluids, as a model for the behavior of "^He.— 
More recently, QMC methods have been used to study 
3D fermionic systems in which there is a hard-sphere re- 
pulsion between opposite-spin particles i^°i^^ To the best 
of our knowledge, however, QMC methods have never 
previously been used to study the phase diagram of 2D 
gases of fermionic hard-core particles 

This paper is arranged as follows. In Sec. |TT] we de- 
scribe the model Hamiltonian that we study. In Sec. 
mil we discuss the form of trial wave function that we 
use. In Sec. IIVI we consider the issue of finite-size errors 
and explain how we extrapolate to the thermodynamic 
limit. In Sec. |V]we present our results for the phase di- 
agram of fermionic hard-core particles and in Sec. IVII we 
present results for the pair-correlation function (PCF), 
static structure factor, and momentum density of hard- 
core fluids. In Sec. IVIII we investigate whether the inclu- 
sion of weak, long-range interactions affects the phase di- 
agram. Finally, we draw our conclusions in Sec. I Villi All 
our QMC calculations were performed using the CASINO 
code.— 



II. HARD-CORE MODEL 

A. Hard-core Hamiltonian 

Suppose we have a system of interacting, circular hard- 
core particles moving in two dimensions. Throughout, we 
use units (a.u.) in which the Dirac constant, the mass of 
the particles, and the radius of the circle that contains 
one particle on average are unity. We assume the par- 
ticles to be fermions with spin s = 1/2. Let D be the 
diameter of the particles. In these units the Hamiltonian 
is 



where 



... if r > D 
I oo otherwise 



(1) 



(2) 



When we refer to "high density," we mean that the value 
of D is large (comparable with the radius of the circle 
that contains one particle on average). In our calcula- 
tions we used finite numbers N of particles in cells subject 
to periodic boundary conditions. From these we extrap- 
olated to infinite system size, as discussed in Sec. IIVI 

The potential energy is zero throughout the permitted 
region of configuration space. Hence the expected po- 
tential energy is zero, and the system only has kinetic 
energy. Nevertheless, the hard-core interaction has a sig- 
nificant effect on the energy, since it defines the boundary 
conditions on the solution of the Schrodinger equation. 
Hence the system does not behave like a noninteracting 
Fermi gas: the magnetic behavior is nontrivial and the 
system must crystallize at a sufficiently high density. 

Hard-core systems have a maximum density: the close- 
packing limit. In 2D the triangular lattice has the highest 
packing fraction, with the maximum hard-core diameter 

being D^ax = y^27r/\/3. 



B. Qualitative features of the phase diagram 

For infinitesimal _D, the system resembles a noninter- 
acting Fermi gas. The momentum density is a unit step 
function. The paramagnetic fluid phase is favored, be- 
cause it has the lowest kinetic energy. The crystal is not 
even stable as an excited state for small D. 

As D is increased, the energy of the paramagnetic 
fluid rises more rapidly than that of the ferromagnetic 
fluid, because wave-function antisymmetry already keeps 
parallel-spin particles apart. The hard-core interactions 
have a greater effect on the distribution of antiparallel- 
spin particles, with the short-range PCF being forced to 
go to zero, as shown in Sec. IVI Al The momentum den- 
sity of a hard-core fluid with _D > is not a unit step 
function, but it retains a discontinuity at the Fermi wave 
vector kp. The crystal becomes stable as D is increased, 
but with a higher energy than that of the fluid. 

At large D, the energies of the ferromagnetic and (frus- 
trated) antiferromagnetic crystals are very similar. At 
some value of D the ferromagnetic fluid is lower in energy 
than the paramagnetic fluid, but at some other value of 
D the crystal becomes more stable than either of the fluid 
phases. It is a nontrivial problem to determine which of 
these transitions occurs first. The crystal orbitals become 
delta functions in the close-packing limit. 



III. TRIAL WAVE FUNCTIONS 

A. Slater-Jastrow-backflow wave functions 

We used trial wave functions of Slater- Jastrow- 
backfiow form ^I* = exp( J)S'^S'-'-. The Jastrow expo- 
nent J included polynomial and plane-wave functions of 
the interparticle distances^^ together with a two-body 
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term to impose the hard-core boundary conditions, as 
described in Sec. IIIIBi and a three-body termi^ The 
backflow function consisted of a polynomial expansion in 
the interparticle distance<^ The terms in the Jastrow ex- 
ponent and backflow function arising from parallel-spin 
and antiparallel-spin pairs of particles were allowed to 
differ. 

Plane- wave orbitals exp(ik • r) were used in the Slater 
determinants S"^ and for the fluid phases. The 
fluid phases suffer from momentum quantization effects 
(single-particle finite-size errors). To reduce these, we 
performed twist averaging in the canonical ensemblci^ 

We used Gaussian orbitals exp (— C|r — RpP) centered 
on hexagonal lattice sites {Rp}, where the exponent C in 
the crystal orbitals is an optimizable parameter, for the 
crystal phases. 

These calculations are similar to the QMC calcula- 
tions that have been performed to establish the zero- 
temperature phase diagram of the three dimensional ho- 
mogeneous electron gas {REG)^^22^ and 2D HEG^"^ 



B. Hard-core two-body behavior 

1. Antiparallel spins 

Let us rewrite the Schrodingcr equation for two hard- 
core particles in terms of the center-of-mass and differ- 
ence coordinates. We assume the center of mass is in its 
zero-energy ground state. The Schrodinger equation for 
the difference coordinate r is 



(3) 



For the boundary conditions we assume that ^' = at 
r = D and — at r ~ R, where R is very large. 
By increasing R, the ground-state eigenvalue E can be 
made arbitrarily close to 0. 5* is circularly symmetric in 
the ground state for antiparallel-spin particles, so 



1 d-^ 
r or 



(4) 



with general solution 5* « A + B\og{r). This approxima- 
tion is valid over any range of r that is small compared 
with R — D. However, we are interested in particle sepa- 
rations of the order D and slightly larger. Applying the 
boundary condition '^{D) = gives 



2. Parallel spins 

Suppose the two hard-core particles have parallel spins. 
In the lowest-energy state the difference-coordinate wave 
function is of the form = W{r)r cos{9) and the en- 
ergy eigenvalue is zero in the limit that the region over 
which the wave function is normalized is large. Hence 
the Schrodingcr equation for the radial part is 



d^W 3 dW 



r dr 



dr"^ 



0. 



(6) 



The general solution over a range of r values that is small 
compared with the region over which the wave function is 
normalized is W(r) « A + Br^^. Applying the boundary 
condition W{D) = 0, one obtains W{r) oc 1 — _D^/r^. So, 
for small r. 



* oc 



1 - ^ ) r cos(6'). 



(7) 



suggesting the parallel-spin two-body Jastrow factor in a 
hard-core gas should be 1 — B"^ jr"^ . 



C. Hard-core Jastrow factor 

The wave function of a hard-core system must go lin- 
early to zero as the separation of any pair of particles 
approaches B. To impose this behavior, the following 
term was included between all pairs of particles in the 
Jastrow exponent J (in addition to the polynomial and 
plane- wave terms): 



UR{r) = \ log 




tanh 



rjP-l 
a(l-r/Lws) 



\ir <B 

\i B <r < Lws 

ff r > iws 



^{r) (X \og{r/B), 



(5) 



(8) 

where L-yvs is the radius of the circle inscribed in the 
Wigncr-Seitz cell of the simulation cell. The parameter 
a was fixed at 1 . uh goes smoothly to zero at Lvvs • The 
other terms in the Jastrow factor are analytic at r = D. 

cxp{uh) ~ taiih.[{r — B) / B] is not of the formlog(r/D) 
or 1 — B^/r^ suggested by the analytic results for 
opposite-spin or same-spin hard-core particles, respec- 
tively, but the short-range behavior is correct [linear in 
(r — B) in the vicinity of r = Z?]. When optimized, the 
polynomial and plane-wave terms in the Jastrow expo- 
nent describe all two-body correlations. Optimized two- 
body Jastrow factors are plotted in Fig. [TJ confirming 
that the simple theory of two-body correlations described 
in Sec. lHIBl is approximately valid. 



suggesting the antiparallel-spin two-body Jastrow factor 
in a hard-core gas should be log(r/D). 

The two-body behavior of hard-core gases was studied 
in Ref. |32| . but this work did not give the two-body Jas- 
trow factor for the 2D gas and did not consider parallel- 
spin pairs. 



D. Free-particle limit 

It is interesting to consider the approach to the free- 
particle limit as the diameter B tends to zero. In a three- 
dimensional hard-core gas, similar arguments to those 
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FIG. 1: (Color online) Two-body Jastrow factors exp[?i(r)] 
for = 50 paramagnetic hard-core particles with D — 0.5 
a.u. 



given m Sees. IIII B il and IIII B 2] show that the two-body 
Jastrow factors are approximately 1~ D/r and 1 — D'^ /r^ 
for antiparallel and parallel spins, respectively, both of 
which tend to unity as D — 0. Hence a Slater- Jastrow 
wave function for a three-dimensional hard-core fluid re- 
duces to a Slater determinant of plane-wave orbitals in 
the limit of zero particle diameter. In the 2D hard- 
core gas the two-body Jastrow factors are approximately 
log(r/_D) and l — D'^/r^ for antiparallel and parallel spins, 
respectively. The antiparallel-spin two-body Jastrow fac- 
tor is the marginal case in which two-body correlations 
become negligible over any given length scale as D is 
made small so that free-particle behavior is recovered in 
the Z? — > limit. In one dimension, however, the in- 
finite contact potential that remains when the _D — > 
limit is taken prevents the hard-core gas from exhibiting 
free-particle behavior, unless the system is fully ferro- 
magnetic. 



coordinate is 
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(9) 



The local energy diverges when hard-core particles ap- 
proach one another, and the sign of the divergence de- 
pends on the positions of all the particles. Hence one 
cannot remove the divergence using a two-body Jastrow 
factor. By contrast, when two charged point particles 
approach one another, the divergence in the local energy 
can be removed by imposing the Kato cusp conditions 
via the two-body Jastrow factor;^ 

The divergence is in principle no worse than that which 
occurs at any other node. However, there are many extra 
nodes introduced by the hard-core potentials, on which 
the local kinetic energy diverges. These nodes can cause 
DMC population explosionsj^i especially if three-body 
terms are omitted from the Jastrow factor. 



F. Need for a three-body Jastrow term 

Physically, it is obvious that at high densities multi- 
body correlation effects will be important: motion is only 
possible if the particles move collectively. In fact, as 
shown in Table HI three-body Jastrow terms lower the 
VMC energy more than backflow. This differs from the 
behavior found in the HEG<^ We have therefore used 
three-body terms in all our calculations. If the range 
of the three-body terms is restricted then the wave func- 
tion becomes significantly poorer, as can be seen in Table 
m Unfortunately, the need to include long-range three- 
body terms in the Jastrow factor makes the calculations 
expensive. 



G. Spin-dependence in paramagnetic phases 



E. Behavior of the local energy as hard-core 
particles collide 



Suppose two antiparallel-spin hard-core particles 1 and 
2 approach each other, i.e., their separation approaches 
D. Their contribution to the Jastrow exponent is J12 = 
log(ri2 — D). Let us write the trial wave function as 
^ = exp( Ji2)S', where S is well-behaved as the particles 
approach. Suppose that all coordinates in the system 
are frozen, apart from the separation of particles 1 and 
2. The contribution to the local energy arising from this 



We used separate two-body Jastrow and two-body 
backflow terms for parallel-spin and antiparallel-spin 
pairs of particles in all our calculations. We considered 
two possible spin-dependences for the three-body terms: 
either (A) using the same three-body term for all triples 
of particles or (B) using separate three-body terms for (i) 
triples involving three particles of the same spin and (ii) 
triples involving two particles of one spin and one par- 
ticle of the opposite spin. VMC results obtained using 
these two different possibilities are shown in Table HIl It 
is clear that allowing different three-body terms for dif- 
ferent spin configurations lowers the variational energy 
significantly (although not by nearly as much as lifting 
the restriction on the range of the three-body terms) and 
hence we have used spin-dependence B in all our produc- 
tion calculations. 
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Meth. 


WF 


SR BF SR 3BJ E (a.u./part.) 


Var. (a.u.) 


VMC 


SJ 






4.8242(4) 


286 


VMC 


SJB 


F 




4.7172(9) 


242 


VMC 


SJB3 


T 


T 


4.6360(2) 


160 


VMC 


SJB3 


F 


T 


4.6272(1) 


158 


VMC 


SJ3 


- 


F 


4.58672(8) 


89 


VMC 


SJB3 


T 


F 


4.5201(2) 


70 


VMC 


SJB3 


F 


F 


4.5072(6) 


65 


DMC 


SJ 


F 


F 


4.4157(6) 




DMC 


SJB 


F 




4.4058(7) 




DMC 


SJB3 


F 


F 


4.3930(8) 





TABLE I: Non-twist-averaged VMC and DMC results with 
different wave functions ("WF") for a paramagnetic fluid of 
N = 26 hard-core particles of diameter D = 1 a.u. The 
"SR BF" column specifies whether or not the backflow func- 
tion was restricted to be short-ranged (cutoff length 2.5 a.u.), 
while the "SR 3BJ" column specifies whether or not the three- 
body Jastrow terms were restricted to be short-ranged (cut- 
off length 2.5 a.u.). Where the cutoff lengths were not re- 
stricted, they approached the radius of the circle inscribed in 
the Wigner-Seitz cell of the simulation cell. "T" and "F" de- 
note true and false, while "Var." is the variance of the local 
energy. 



Spin- 


-dep. SR BF & 3BJ Energy (a.u./part.) 


Var. (a.u.) 


A 


T 


4.640(1) 


165 


B 


T 


4.6360(2) 


160 


A 


F 


4.534(1) 


76 


B 


F 


4.5072(6) 


65 



TABLE II: Non-twist-averaged VMC results with different 
spin dependences for a paramagnetic fluid of A*" = 26 hard- 
core particles of diameter D = 1 a.u. The "SR BF & 3BJ" 
column specifies whether or not the backflow function and 
three-body Jastrow terms were restricted to be short-ranged 
(cutoff length 2.5 a.u.). Where the cutoff lengths were not 
restricted, they approached the radius of the circle inscribed 
in the Wigner-Seitz cell of the simulation cell. "T" and "F" 
denote true and false, while "Var." is the variance of the local 
energy. 



H. Relative accuracy of wave functions for 
different phases 



The variance of the energy is zero if the trial wave func- 
tion is an eigenfunetion of the Hamiltonian. As shown in 
Table IIIIl the variance per particle is significantly lower 
for the ferromagnetic fluid than the paramagnetic fiuid, 
indicating that the trial wave function is more accurate 
for the former than the latter. The variance per parti- 
cle is similar for ferromagnetic fluids and crystals near 
the transition density. Hence, if anything, our results are 
biased in favor of ferromagnetic phases. 



D (a.u.) 


1 Phase 


Var. per part, (a.u.) 


1 


Para, fluid 


2.49 


1 


Ferro. fluid 


1.76 


1 


Ferro. crystal 


0.736 


0.88 


Para, fluid 


1.34 


0.88 


Ferro. fluid 


0.616 


0.88 


Ferro. crystal 


0.494 


0.63 


Para, fluid 


0.385 


0.63 


Ferro. fluid 


0.147 


0.63 


Ferro. crystal 


0.228 



TABLE III: Non-twist-averaged VMC results for paramag- 
netic fluids with A'' = 26 hard-core particles and ferromag- 
netic fluids with A = 25 hard-core particles. 



I. Nature of phase transition 

We have looked for a first-order phase transition by 
comparing fixed-node DMC energies with fluid and crys- 
tal orbitals, relying on the flxed-node approximation to 
impose the symmetry of the phase on the wave func- 
tion. However, it is possible that there could actually 
be a continuous transition from fluid to crystal behavior. 
The fact that the "fluid" wave function tries to become 
crystal-like at high density and the "crystal" wave func- 
tion tries to become fluid-like at low density (see Sec. lVIj) 
supports this view. Nevertheless, even if this is the case, 
our calculations determine the region in which crystal- 
lization is expected to take place, and demonstrate that 
a ferromagnetic fluid phase is unlikely to occur. 



IV. FINITE-SIZE EFFECTS 



According to the simple theory given in Sec. IHI B 2\ 
the long-range parallel-spin two-body Jastrow exponent 
is approximately given by 

Ucair) = log(l - D^/r^) « -D^/,-^ + 0{r-^). (10) 

The 2D Fourier transform of the leading term does not 
have a power-law behavior. In fact numerical tests show 
that limjj^o / exp(zk • r — ar) dr diverges logarithmi- 
cally at small k. So, for small fc, Uaa{k) may be written 
as —cD^ log(fc), where c is a positive constant. 

The Chiesa-Holzmann-Martin-Ceperle}*^^ approxima- 
tion to the long-range two-body finite-size correction to 
the kinetic energy in 2D is 



AT 



E 



4(27r)2 



2'Kk X k'^Uaa{k) dk 



= / k\^^{k)dk, (11) 

where Q — 2^/ttJA is the radius of the circle in k-space 
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with area (27r)^/A and A = ttN is the area of the simu- 
lation ceU. Inserting Uaa{k) = — cD^ log(fc) gives 



AT 



ttNcD^ 
0[D^N-^ \ogiN)] 



log (2^yV/Tj - Att^/A^ 



(12) 



Hence, ignoring the logarithmic factor, the leading-order 
correction to the energy per particle due to the neglect 
of long-range two-body correlations falls off as 0{N~'^). 
This is much more rapid than the analogous result for 
the 2D HEG.— The leading-order correction is positive. 
Unsurprisingly, however, two-body finite-size errors get 
rapidly more severe as D increases. 

We have carried out numerical tests which show that 
the bias due to twist averaging in the canonical en- 
semble in two dimensions falls off as 0{N~'^/'^), but 
with a small prefactor. Residual single-particle errors in 
the canonical-ensemble twist-averaged fluid energies are 
small compared with the effect of twist averaging, and 
the difference between the twist-averaged and non-twist- 
averaged energies is only about 0.02 a.u. per particle at 
D = l a.u. for iV = 10 and 26. Hence the 0{N~'^) behav- 
ior due to the neglect of long-range two-body correlations 
dominates the systematic error for all the hard-core di- 
ameters D that we have studied, as can be seen in Fig. 

m 

Several data points in Fig. [5] are outliers. This non- 
systematic behavior appears to be a genuine finite-size 
effect. The finite-size "noise" in the crystal energies is 
more severe at low densities, while the finite-size noise 
in the fluid energies is more severe at very high densi- 
ties. The theory of two-body finite-size effects developed 
above breaks down for fluids at high density and crystals 
at low density: the sign of the bias is wrong. These cases 
are pathological in various respects. At very high den- 
sities, the fluid energies obtained with the same number 
of particles in different-shaped cells disagree, although 
they extrapolate to the same value in the limit of infinite 
system size. In fact at very high density it is geomet- 
rically impossible to fit finite fluids in some cell shapes. 
As can be seen in Fig. [3l the Gaussian exponent C of the 
crystal is not well-behaved for D < 0.75, suggesting that 
the crystal is becoming unstable. As shown in Fig. [T4l 
the momentum density is developing an edge, leading to 
substantial single-particle finite-size errors. 

As shown in Table U restricting the range of the three- 
body term raises the variational energy significantly. 
Three-body correlations arc therforc long-ranged. So 
there must be three-body finite-size errors in the fluid 
and crystal energies obtained in finite simulation cells. 
We assume they are a "random" error about the system- 
atic finite-size bias due to two-body finite-size effects. We 
have therefore obtained several E{N) data points for each 
phase at each density in order to average out the "noise" 
when extrapolating to inflnite system size. 
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FIG. 2: (Color online) DMC energy E against the reciprocal 
of the square of the system size A''"'^ for (a) D = 0.5 a.u., (b) 
D = 0.63 a.u., (c) D = 0.75 a.u., (d) D = 0.88 a.u., and (e) 
D = 1 a.u. 



V. PHASE DIAGRAM 

The energies of the different phases are compared in 
Fig. U and are plotted relative to the energy of the para- 
magnetic fluid phase in Fig. [51 The continuous curves 
shown arc Akima spline interpolations between the DMC 
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FIG. 3: (Color online) Optimal Gaussian exponent C for the 
ferromagnetic crystal phase against the reciprocal of system 
size iV"^. 
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FIG. 5: (Color online) DMC energy of the different phases 
relative to the energy of the paramagnetic fluid. 
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FIG. 4: (Color online) DMC energy E against hard-core di- 
ameter D. Note that, when D = the energies of the param- 
agnetic and ferromagnetic fluids are 0.5 and 1 a.u. per particle, 
respectively. 



VI. OTHER PROPERTIES OF HARD-CORE 
GASES 

A. Pair-correlation function 

We have calculated the PCFs of the fluid phases of the 
hard-core gas, and our results arc shown in Fig. |6l The 
PCFs, which were obtained without twist averaging, are 
not well-converged with respect to system size at high 
density: the finite-size errors are oscillatory. Neverthe- 
less, we can make some qualitative comments about the 
physics revealed by the PCF. 

The distance between the peaks of the PCF obtained 
in square cells at 13 = 1 a.u. is y/ir, which arises from the 
square-cell geometry. The difference between the PCFs 
in hexagonal and square cells is still significant at 13 = 
0.88 a.u. and = 26. The antiparallel-spin PCF evolves 
slowly towards g{T) = 1 (the free-particle result) as D is 
reduced. The size of the exchange-correlation hole and 
the ripples in the parallel-spin PCF are determined by 
the Fermi wave vector kp when D is small. 



energies. (Akima spline interpolation is stable to the 
presence of outliers in the datai^) It can be seen that 
crystallization takes place when D = 0.83 a.u., leaving 
no region of stability for a ferromagnetic fluid. (Recall 
that our calculations are, if anything, biased in favor of 
the ferromagnetic fluid.) Thus our calculations rule out 
the possibility of an itinerant ferromagnetic fluid phase 
in a 2D gas of particles with only hard-core interactions. 
For those values of D for which the crystal has the lowest 
energy, the energy difference between the antiferromag- 
netic and ferromagnetic states is insignificant, showing 
that exchange interactions in the crystalline phase are 
negligible. 



B. Static structure factor 

The static structure factor of a paramagnetic fiuid at 
D — 0.63 a.u. for various values of N is plotted in Fig. 
[T] The structure factor is well converged with respect to 
N. 

VMC-calculated structure factors for paramagnetic 
and ferromagnetic hard-core gases are shown in Figs. [8] 
and [HI respectively. Finite-size "noise" is much worse 
for ferromagnetic phases. The fiuid and crystal structure 
factors arc similar, especially at low densities. 
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FIG. 6: (Color online) VMC PCFs g{r) of hard-core systems 
at (a) D = 0.5 a.u., (b) D = 0.75 a.u., (c) D = 0.88 a.u., and 
(d) D — 1 a.u., for different system sizes A^. 
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FIG. 7: (Color onhne) VMC structure factor of the param- 
agnetic fluid at D = 0.63 a.u. for a range of system sizes A*'. 




FIG. 8: (Color online) VMC structure factors of paramagnetic 
fluids. For D — the analytic result for a paramagnetic free- 
particle gas is given. 
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FIG. 9: (Color online) VMC structure factors of ferromag- 
netic fluids and crystals. The D = result is the analytic 
structure factor for a ferromagnetic free-particle gas. 
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FIG. 10: (Color online) VMC momentum density of the para- 
magnetic fluid at D = 0.88 a.u. at different system sizes A^. 



FIG. 12: (Color online) VMC momentum densities of ferro- 
magnetic fluids at different hard-core diameters D. 
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FIG. 11: (Color online) VMC momentum densities of param- 
agnetic fluids at different hard-core diameters D. 
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FIG. 13: (Color online) Renormalization factor Z against D, 
from fits to the momentum density. 



C. Momentum density 

1. Results for fluids 

The momentum density shown in Fig.[TU]is reasonably 
well converged with respect to N. VMC results for the 
momentum densities of paramagnetic and ferromagnetic 
hard-core systems are shown in Figs. [TT] and [T^l respec- 
tively. 



2. Renormalization factors 

The renormalization factor — the discontinuity Z in the 
momentum density at the Fermi edge — is plotted against 
hard-core diameter in Fig. 1131 The finite-size errors in 
Z arc oscillatory, so we have averaged over the results 
obtained for different system sizes N. 



3. Results for crystals 

The momentum densities of ferromagnetic hard-core 
crystals arc shown in Fig. [TH This figure strongly 
confirms the conclusion that the crystal is unstable at 
D = 0.63 a.u.: the "crystal" momentum density develops 
a near-discontinuity at the Fermi wave vector. 



VII. WEAK, LONG-RANGE INTERACTIONS 

A. Form of the interaction 

Our results show that, close to the transition from fluid 
to crystal, the energies of the paramagnetic and ferro- 
magnetic fluids are very finely balanced, albeit with the 
paramagnetic fluid always lying lower in energy. It is in- 
teresting to ask whether a small change in the two-body 
potential might alter the relative stabilities of these two 
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FIG. 14: (Color online) VMC momentum densities of ferro- 
magnetic crystals. 



FIG. 15: (Color online) Change in energy resulting from tiie 
inclusion of a Ar~^ interaction for a paramagnetic fluid of 
= 18 hard-core particles of diameter D = 0.88 a.u. 



fluid phases and allow a region of stable ferromagnetic 
fluid. 

As explained in Sec. U it can be arranged that 
fermionie molecules confined to a plane experience an 
additional long-range, pairwise interaction varying as 
Ar~^, where r is the interparticle separation and A is a 
constantj^ii In a finite, periodic cell, the two-body inter- 
action between one particle and all the images of another 
particle that is at a distance r from the first is 
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3f2 



4i?| 



(13) 



where S" is a circular region of radius Rs centered on the 
origin, and f is the minimum image of r. By making Rs 
sufficiently large, the approximation to the infinite sum 
can be made arbitrarily good. 

The "Madelung" term (the interaction of each particle 
with its own images) is 



Rt^O 



A 



E 



RGS-{0} 



A_ 

m 



A 



E |R|3 

ReS-{o} ' ' 



R>Rs 

27rA 



AR 



s 



(14) 



The Hamiltonian for the finite, periodic cell is therefore 
H = -\y^l + Y.M\hA) + VA{h,)] + ^. (15) 

i i>j 



B. Convergence of the real-space sum 

We have chosen Rg such that 119 stars of lattice vec- 
tors are summed over explicitly in square cells (for fluid 
phases) and 86 stars are summed over in hexagonal cells 
(for crystal phases). The error in the Madelung constant 
is about 2 x 10^*" a.u. per particle for iV = 18 particles 
(and smaller at larger N). The error in the additional 
interaction per particle due to the finite number of stars 
of R vectors should therefore be much smaller than 10~^ 
a.u. The error made in truncating the real-space sum is 
therefore negligibly small compared with the statistical 
error bars on our QMC energies. 



C. Validity of perturbation theory 

To a very good approximation we can describe the ef- 
fect of the Ar~3 interaction using perturbation theory. 
The change in the energy resulting from the additional in- 
teraction is approximately given by the expectation of the 
interaction operator with respect to the previously opti- 
mized wave function for the unperturbed system, which 
can be evaluated using VMC. For — 1 < A < 1, this ap- 
proximation reproduces the full DMC energy difference 
to within about 0.002 a.u. at = 0.88 a.u., as shown in 



D. Finite-size errors 

As with total energies, finite-size errors in the expec- 
tation value of the extra interaction are larger and quasi- 
random for the unstable phases. We extrapolate to infi- 
nite system size in each case, as shown in Fig. 1161 
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FIG. 16: (Color online) VMC-evaluated expectation value 
of ^iyj r^j^ with respect to the non-reoptimized VMC wave 
function. 



E. Effect of weak interactions on the phase diagram 



At D = 0.88 a.u., combining the VMC results for 



with the DMC energies in the absence of the 



additional interaction, we find that 

= 3.2559 + 0.9296A 
Eff = 3.251 + 0.928A 
-Bfc = 3.2305 + 0.9244A, 



(16) 



where iJpF, i^PF, and E-pc are the energies in a.u. per par- 
ticle for the paramagnetic fluid, ferromagnetic fluid, and 
ferromagnetic crystals, respectively, as a function of A. 
The resulting offset in the energy relative to the energy 
of the paramagnetic fluid phase at a density close to the 
crystallization density is shown in Fig. [171 The changes 
in the energies of the three phases due to the inclusion 
of the A.r~^ tail are very similar. As one would expect, a 
repulsive potential (A > 0) favors phases where the parti- 
cles are kept apart (i.e., ferromagnetic over paramagnetic 
phases, and crystals over fluids), whereas an attractive 
interaction (A < 0) favors a paramagnetic fluid the most 
and a ferromagnetic crystal the least. There is no region 
of stability for the ferromagnetic fluid, although it comes 
close at A k, —5. However, this is a regime where the 
interaction is strongly attractive, so that perturbation 
theory is no longer valid, and a trial wave function that 
includes the possibility of superfluid pairing is required 
for accurate QMC calculations. It therefore seems un- 
likely that a weak interaction will have much effect on 



VIII. CONCLUSIONS 



Gases of hard-core particles are a natural model for 
ultracold-atom systems. Correlation effects in 2D hard- 
core systen is are surprisingly long-ranged. Our Q MC 
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FIG. 17: (Color online) Energy of the ferromagnetic fluid and 
crystal relative to the energy of the paramagnetic fluid against 
the strength of the interaction for hard-core particles of 
diameter D = 0.88 a.u. 



results show that there is no regime in which itinerant 
ferromagnetism occurs in a 2D hard-core fluid. As the 
hard-core diameter D is increased, the system undergoes 
a transition from a paramagnetic fluid to a crystal at 
D ~ 0.83 a.u. The absence of a region of stability for 
a ferromagnetic fluid resembles the situation in the 2D 
HEGj^ Including a weak tail in the two-body in- 
teraction between hard-core particles does not lead to a 
significant revision of the phase diagram. We have pre- 
sented QMC results for the PCFs, static structure fac- 
tors, and momentum densities of the fluid phases of the 
hard-core gas. 
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